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I .  INTRODUCTION 


Advances  in  the  development  of  laser  based  combustion  diagnostic 
techniques,  of  numerical  algorithms  for  solving  ordinary  and  partial 
differential  equations,  and  of  larger  capacity  and  higher  speed  computing 
machines  have  combined  to  enable  one  to  begin  to  construct  models  of  com¬ 
plex  combustion  phenomena  that  include  fundamental  physico-chemical  pro¬ 
cesses.  The  necessary  inputs  to  such  models  are  data  sets  of  (1) 
chemical  reactions,  including  their  rate  coefficients,  (2)  thermodynamic 
coefficients  and,  (3)  transport  coefficients. 

In  the  work  presented  here,  we  have  used  a  one-dimensional,  premixed, 
laminar,  steady-state  flame  model  to  predict  individual  species  and 
temperature  profiles.  The  computed  profiles  are  compared  with  experi¬ 
mental  data  to  test  the  model.  Essential  agreement  between  measured 
and  computed  profiles  is  taken  as  validation  of  the  model  input  parameters. 
On  the  other  hand,  discrepancies  indicate  those  areas  that  are  in  need 
of  further  study. 

We  have  studied  the  ozone  flame  as  a  test  case  because  of  its 
simplicity.  There  are  three  species:  0,  02,  and  O3,  related  by  the 
reactions 


0^  +  M  +  0  +  0^  +  M, 
k-i 


(i) 


0  +  °3  t  °2  +  °2» 
k-2 


(II) 


and 


0  +  0  + 


M  t 


CL  +  M. 


-3 


(III) 


We  emphasize  that  the  requisite  input  data  are  obtained  from  independent 
experiments  and  that  adjustable  (arbitrary)  parameters  or  average  values 
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were  not  employed.  Specifically,  the  thermodynamic,  kinetic  and  laminar 
transport  properties  are  both  species  and  temperature  dependent. 

With  few  exceptions,1 2"3 4 * * 7  most  other  models  of  the  ozone  flame  have 
emphasized  mathematical  approaches  or  numerical  techniques . 4-7  That  is, 
their  primary  objective  has  been  the  development  of  methods  of  solution 
rather  than  detailed  comparisons  of  model  predictions  and  experimental 
results.  Consequently  these  authors  either  employed  mathematically  non- 
essential  assumptions  that  simplified  the  transport  or  employed  kinetic 
coefficients  that  are  incompatible  with  what  is  presently  known,  or  both 
They  have  also  generally  limited  their  illustrations  to  the  single  case 
of  0.25  initial  ozone  mole  fraction  (IOMF). 

The  object  of  this  paper  is  to  present  as  detailed  an  assembly  of 
input  parameters  as  practical,  to  describe  the  input  coefficients  by 
means  of  convenient  expressions,  to  compare  our  model  results  with  other 
results  and  to  indicate  future  directions  of  investigation.  We  do  this 
so  that  a  complete  verified  set  of  realistic  input  parameters  may  become 
available  for  what  has  become  a  standard  test  case,  the  premixed  laminar 
ozone  flame. 


F.  Cramarossa  and  G.  Dixon-Lewis ,  "Ozone  Decomposition  in  Relation  to 
the  Problem  of  the  Existence  of  Steady-State  Flames ",  Combustion  and 
Flame  16,  243-251  (1971). 

2 

Kenneth  A.  Wilde,  "Boundary -Value  Solutions  of  the  One-Dimensional 
Laminar  Flame  Propagation  Equations",  Combustion  and  Flame  18,  43-52, 
(1972). 

£ 

J.  Warnatz,  "Calculation  of  the  Structure  of  Laminar  Flat  Flames  I: 

Flame  Velocity  of  Freely  Propagating  Ozone  Decomposition  Flames", 

Ber.  Bunsenges  Phys.  Chem.  82,  193-200  (1978). 

4 

Leon  Bledjian,  "Computation  of  Time-Dependent  Laminar  Flame  Structure" , 
Combustion  and  Flame  20,  5-17 ,  (1973). 

Edwin  S.  Campbell,  "A  Theoretical  Analysis  of  Chemical  and  Physical 
Processes  in  an  Ozone  Flame",  Chem.  Engineering  Sci.  20,  311-329,  (1965). 

Q 

Stephen  B.  Margolis ,  "Time-Dependent  Solution  of  a  Premixed  Laminar 
Flame",  J.  of  Computational  Phys.  27,  410-427,  (1978). 

7 

J.  0.  Hirschf elder ,  C.  F.  Curtis  and  Dorothy  E.  Campbell,  "The  Theory 
of  Flame  Propagation.  IV",  J.  Phys.  Chem.  57,  403-414,  (1953). 
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II.  THE  FLAME  EQUATIONS 

The  derivation  of  the  conservation  equations  for  a  multi -component 
reacting  ideal  gas  mixture  can  be  found  in  the  literature. 8-10  Here  we 
are  interested  in  those  equations  that  adequately  describe  a  one¬ 
dimensional,  laminar,  premixed  flame  that  propogates  in  an  unbounded 
medium.  The  effects  of  radiation,  viscosity,  thermal  diffusion  and  body 
forces  are  ignored.  Since  the  burning  velocity  is  small  compared  with 
the  local  speed  of  sound,  the  pressure  is  taken  to  be  constant . 10, 11 

The  equations  pertinent  to  our  present  case  are: 

Overall  Continuity: 


ie.  +  Hpu)  _  0 
at  3x  "  u 


Ci) 


Continuity  of  Species: 

3Yk  3Yk  3 

P  3^*  pu  53T=  ■  lx  tpYkV  *  RkMk>  k  -  1.2.... N  (2) 

and 

Conservation  of  Energy: 


pc  9T  puc  3T 


P  8t 


3x 


3_ 

3x 


N 

E 

k=l 


N  3T 

R,  M,  h.  -p  E  c  ,  Y  V  — 
k  k  k  j.  pk  k  k  3x, 


(3) 


where  the  variables  are  defined  in  the  glossary.  In  equations  (2)  and 
(3),  Vk  is  determined  by  the  Stefan-Maxwell  equations. 


Q 

J.  0.  Hirschf elder,  C.  F.  Curtis  and  R.  B.  Bird,  Molecular  Theory  of 
Gases  and  Liquidsj  2nd  Printing,  Corrected,  with  notes,  John  Wiley  and 
Sons,  NY,  (1964).  u 

9 

Bird,  W.  S.  Stewart  and  E.  N.  Light  foot.  Transport  Phenomena , 
John  Wiley  and  Sons,  NY,  (1960).  '  f 

°F.  A.  Williams,  Combustion  Theory  Addison-Wesley ,  Reading,  MA,  (1965). 

R.  M.  Fristrom  and  A.  A.  Westenberg,  Flame  Structure ,  McGrccw-Hill  NY 
(1965)  p.  319.  - ** 
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(4) 


N 

SA  Dkj 


X,  X. 

Ju. 


The  boundary  conditions  are  the  following  (t  >  0) . 

For  x  =  +  00 

T  =  T  and 
u 

Yk  ■  Yku-  <k  - 
and  for  x  =  -  00 


9T  _  _ k  _  ,,  _  .  -  ,n 

3x  3x  °»  (k  ~  1>2> • •  -N) . 


In  order  to  avoid  solving  equation  (1)  a  coordinate  ip  is  introduced 
such  that 


x 

<Kx»t)  =  /  p(x'  ,t)  dx'. 
o 


(6) 


Then  If  =  p  and  If  =  -pu  +  mo  Where  pul  x=0  =  mo 

With  this  notation  equations  (2)  and  (3)  become 


3Y.  3Y. 

k  .  k 

-  +  m  - 

at  o  dip 


fir  c»W 


RkMk/p 


and 


(7) 


3T 

3t 


+  m 


3T 

O  dip 


1_ 

C 

P 


3  ,  .  3T. 

3^  (-pX34/-) 


1_ 

pc 


N 

E 

k=l 


RkMkhk 


N 

E  !pk  PY  V  dT  (8) 

k=l  C  K  dip 

p 
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respectively.  The  solution  of  equations  (7)  and  (8)  yields  the  steady 
state  profiles  of  atomic  oxygen,  ozone  and  temperature.  (In  practice 
convenient  dimensionless  forms  of  equations  (7)  and  C8)  are  solved.) 
The  molecular  oxygen  profile  is  obtained  from  the  relation 


3 

I 

k=l 


1 . 


The  burning  velocity,  Su,  is  defined  as  the  velocity  of  the  flame 
relative  to  the  fluid  at  rest,  i.e.,  at  infinity.  The  value  of  Su  can 
be  found  from  the  steady  state  profiles  by  performing  a  coordinate 
transformation  such  that  the  flame  appears  stationary.  In  this  new 
system  all  variables  are  independent  of  time;  specifically. 


8Yk 

=  0,  and 

=  0. 

From  equation  (1)  we  then  have 

8(pu)/8x  =  0,  or  pu  =  constant. 

Take  any  one  of  the  equations  (2)  and  integrate  over  any  interval  (a,b) 
to  obtain 


8t 

9p 

at 


pu[Yk(b)  -  Yk(a)] 


b  b 

!  Wx  -  pYkVk| 


Then 


u  (°°) 


b  b 

'a  Wx  -  eYkVa 

p(“)[Yk(b)  -  Yk(a)] 


b  b 

U  P-1  W*  '  pYkVkla  . 

P(“HYk(b)  -  Yk(a)] 


In  this  coordinate  system  Su  =  -  u  (°°)  . 


(9) 


(10) 
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III.  PDECOL  -  A  GENERAL  PDE  SOLVER 

The  package  PDECOL,  developed  by  Madsen  and  Sincovec-^^  was  used  to 
solve  this  problem.  This  is  a  general  package  for  solving  partial 
differential  equations  (PDE's)  using  the  method  of  lines. 

The  spatial  discretization  is  accomplished  by  finite  element  col¬ 
location  methods  based  on  B-splines. 13  The  basic  assumption  is  that  the 
solution  can  be  written  in  the  form 


Yk  =  gC  c^)(t)Bi(i|»),  k=l...N,  and 
i=i 


i=i 


(t)B.  oio. 


where  the  basis  functions  Bj(ip),  i  =  l  , ..NC,  are  B-splines  and  span  the 
solution  space  for  any  fixed  t  to  within  a  small  error  tolerance.  The 
time  dependent  coefficients  are  determined  uniquely  by  requiring 

that  the  expansion  above  satisfy  the  given  boundary  conditions  and  that 
they  satisfy  equations  (7)  and  (8)  exactly  at  (NC-2)  interior  (colloca¬ 
tion)  points.  If  there  is  a  null  boundary  condition,  an  extra  colloca¬ 
tion  point  is  added.  The  Bj  are  piecewise  polynomials  of  order  KORD . 
Given  a  user  supplied  set  of  NB  breakpoints  (i.e.,  a  set  of  strictly 
increasing  locations  where  the  polynomials  are  joined)  and  the  number 
of  continuity  conditions,  NCC,  to  be  applied  at  the  breakpoints,  PDECOL 
generates  a  set  NC  =  KORD  (NB-1)  -  NCC  (NB-2)  basis  functions  and  col¬ 
location  points.  Since  by  definition  a  B-spline  is  zero  except  over  a 
small  interval,  at  any  collocation  point  no  more  than  KORD  of  the  B- 
splines  are  nonzero.  So  the  system  of  0DE!s  for  the  coefficients  CfcCi) 
will  not  be  fully  coupled. 


12 

B.  K.  Madsen  and  R.  F.  Sincovec3  "PDECOL:  General  Collocation  Soft¬ 
ware  for  Partial  Differential  Equations" ,  Preprint  UCRL-78263  (Rev  1) 
Lawrence  Livermore  Laboratory 3  (1977). 

13 

C.  de  Boor3  "Package  for  Calculating  With  B-Splines" 3  Siam.  J.  Numer. 
Anal.  14,  441-472 3  (1977). 
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This  system  of  ODE's  is  integrated  in  time,  using  a  variant  of  the 
Gear  integrator.  The  appropriate  banded  Jacobian  is  generated  internally 
by  the  program.  Once  the  integrator  has  reached  a  desired  output  time  t, 
the  values  of  Yk  and  T  can  be  obtained  for  any  >  by  substituting  into 
the  expansion. 


The  definition  of  the  set  of  breakpoints  that  determines  the  spatial 
discretization  is  crucial.  They  must  be  chosen  close  enough  that  spatial 
errors  do  not  destroy  the  solution  yet  not  so  dense  that  one's  computer 
resources  are  exceeded.  There  is  at  present  no  way  of  determining  an 
optimum  set  of  breakpoints  for  a  given  problem  other  than  trial  and 
error.  We  have  developed!!  a  method  of  concentrating  our  breakpoints  in 
the  steep  flame-front,  where  accuracy  is  necessary.  This  method  is 
similar  to  a  procedure  developed  by  Spalding  and  StephensonlS  for  use  in 
a  finite  difference  code. 


With  initial  guessed  profiles  for  the  ozone  flame  and  with  values 
of  KORD  =  6,  NCC  =  5  and  59  £  NB  <_  70,  the  program  coding  has  executed 
in  30  to  60  seconds  on  the  BRL  Cyber  76  and  has  required  about  80,0008 
words  of  core.  A  listing  of  the  program  is  given  in  reference  14. 


IV.  INPUT  COEFFICIENTS 

A.  Kinetic 

Reaction  (1)  is  the  sine  qua  non  of  the  ozone  flame  and  as  such  is 
the  single  most  important  reaction.  Thus  we  require  the  expression  of 
k}  to  be  valid  over  the  range  of  adiabatic  flame  temperatures,  <v  1000K- 
3000K.  Available  reviewsl6-17  recommend  expressions  that  are  valid  only 
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T.  P.  Coffee  and  J .  M.  Heimerl ,  "A  Method  for  Computing  the  Flame 
Speed  fov  a  Laminar,  Premixed,  One  Dimensional  Flame",  BRL  Technical 
Report ,  ARBRL-TR-0 2212 ,  Jan  80.  (AD#A082803) 

1 5 

D.  B.  Spalding  and  P .  L.  Stephenson,  "Laminar  Flame  Propagation  in 
Hydrogen  and  Bromine  Mixtures ",  Proc.  Roy.  Soc.  London  A3 24,  315-337, 
1971-  (See  also),  D.  B.  Spalding,  P.  L.  Stephenson  and  R.  G.  Taylor, 

"A  Calculation  Procedure  for  the  Prediction  of  Laminar  Flame  Speeds" 
Combustion  and  Flame,  17_,  55-64,  1971. 

1 6 

D\  L.  Baulch,  D.  D..  Rrysdale,  J.  Duxbury  and  S.  J.  Grants  Evaluated 
Kinetic  Data  fov  High  Temperature  Reactions ,  Vo l  3 ,  Homogeneous  gas 
phase  reactions  of  the  O2/O3  system,  the  CO/O2/H2  system  and  of  sulphur 
containing  species.  Butterworths ,  Boston,  (1976). 

17 

Harold  S.  Johnston,  rrGas  Phase  Reaction  Kinetics  of  Heutral  Oxygen 
Specie s!t ,  NSRDS-NBS-20,  September  1968. 
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up  to  1000K.  Recently  we!8  have  critically  examined  the  high  temperature 
shock  tube  results  of  Michael^  (971K-1384K)  and  of  Center  and  Kung20 
(2041K-2941K)  and  found  the  expression  listed  in  Table  1.  This  expression 
describes  all  the  direct  experimental  data  available  for  ki  and  is  valid 
over  the  range  300K-3000K.  To  a  good  approximation  it  is  independent  of 
any  set  of  experimental  data.  (For  a  more  detailed  discussion,  see  re¬ 
ference  18). 

The  expression  for  k2  and  k3  have  been  taken  from  the  review  liter¬ 
ature  of  Hampson2!  and  Johnston, 17  respectively.  (Hampson's  and 
Johnston's  expression  for  k2  yield  values  that  differ  by  15%  or  less 
over  the  range  600K  <_  T  <_  3000K) .  The  expression  of  k-i  and  k_2  have 
been  derived  from  the  forward  rate  coefficients  and  the  corresponding 
equilibrium  constants.  The  expression  for  k-3  is  taken  from  Johnston. 

Over  the  range  1000K  <  T<  3000K  the  value  of  the  equilibrium  constants 
for  reaction  3  as  computed  from  Johnston's  expressions  for  k3  and  k_3 
is  at  most  20%  greater  than  published  equilibrium  values. 16 

The  rate  coefficient  may  vary  with  the  identity  of  the  chaperon  or 
third  body,  M.  The  chaperon  relationships  used  are  ki(M=02)  =  0.44 
kl(M=03),17  k3(M=0)  =3.6  k3(M=02),16  ki(M=0)  =  ki(M=02)  and  k3(M=03)  = 
k3 (M-O2) .  Only  the  first  relationship  has  been  found  to  be  important. 

B.  Thermodynamic 


The  enthalpy  and  heat  capacity  are  functions  of  species  and  of  tem¬ 
perature.  The  JANAF22  values  for  the  allotropes  of  oxygen  are  described 
to  better  than  a  few  parts  per  thousand  over  the  temperature  range  300K- 


1 8 

J.  M.  Heimerl  and  T.  P.  Coffee,  "The  Unimolecular  Ozone  Decomposition 
Reaction ",  Combustion  and  Flame ,  35,  117-123,  1979.  For  a  complete 
listing  of  all  direct  measurements,  see  BRL  Technical  Report  ARBRL-TR- 
02185,  Aug  79.  (AD#A076966) 

19 

J.  V.  Michael,  "Thermal  Decomposition  of  Ozone",  J.  Chem.  Phys.  54, 
4455-4459,  (1971).  — 

20 

R.  E.  Center  and  R.  T.  V.  Kung,  "Shock  Tube  Study  of  the  Thermal  De¬ 
composition  of  O3  From  1000  to  3000°K ",  J.  Chem.  Phys.  62,  801-807 . 
(1975).  — 

21 

R.'  F.  Hampson,  (ed) ,  "Survey  of  Photochemical  and  Rate  Data  for  Twenty- 
Eight  Reactions  of  Interest  in  Atmospheric  Chemistry ",  J.  Phys.  Chem . 
Ref,  Data  2,  267-312,  (1973). 

22 

D.  R.  Stull  and  H.  Prophet,  JANAF  Thermochemical  Tables ,  2nd  Edition. 
NSRDS-NBS-37,  June  1971. 
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TABLE  1.  KINETIC  COEFFICIENTS 
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3000K  by  the  polynomial  fits  of  Gordon  and  McBride.23  These  fits  take 
the  form: 


H°  =  R(a  +  E  a  Tn/n) 
I  o  _  n  J 

n=l 

and 


C  =  R  l  a  Tn_1, 
P  n  n 

1  TT  —  1 


where  the  coefficients  are  listed  in  Table  2. 

C .  Transport 

The  transport  coefficients  used  for  the  ozone  flame  are  listed  in 
lable  3.  (The  subscripts  1,  2,  and  3  for  the  transport  coefficents  stand 
for  0,  02  and  O3  respectively).  Expressions  for  Al,  A2  and  pDi2  have 
been  taken  from  the  literature.  Dalgarno  and  Smith* 24  have  computed  the 
thermal  conductivity  for  atomic  oxygen  up  to  100,000K  and  derive  an 
analytic  expression  for  Al  valid  over  the  range  100K-2000K.  Hanley  and 
Ely^-3  have  critically  reviewed  the  thermal  conductivity  data  for  mole¬ 
cular  oxygen  and  have  tabulated  values  for  A2  over  the  range  80K-2000K. 
Using  the  functional  form  AT^  we  have  obtained  a  fit  to  their  tabular 
values  over  the  range  300K-2000K  which  fit  has  an  overall  root  -mean-square 
error  of  ,vl%  and  a  maximum  error  of  -v2%.  Marrero  and  Mason2^  have 
critically  reviewed  data  for  pDi2  and  determined  an  expression  valid 
over  the  temperature  range  280K-10, 000K.  It  has  an  estimated  uncertainty 
of  ±  25%  for  temperatures  above  1000K. 


S.  Gordon  and  B.  J.  McBride,  "Computer  Program  for  Calculation  of  Com¬ 
plex  Chemical  Equilibrium  Compositions ,  Rocket  Performance,  Incident 
and  Reflected  Shocks,  and  Chapman- J ouguet  Detonations" ,  NASA-SP-273, 
(1971),  (1976  program  version) . 

24 

A.  Dalgarno  and  F.  J.  Smith,  "The  Thermal  Conductivity  and  Viscosity 
of  Atomic  Oxygen ",  Planet  Spac.  Sci.  9^  1-2,  (1962). 

25 

H.  J.  M.'  Hanley  and  J.  F.  Ely,  "The  Viscosity  and  Thermal  Conductivity 
Coefficients  of  Dilute  nitrogen  and  Oxygen",  J.  Phys.  Chem.  Ref.  Data, 

2,  735-755,  (1973). 

2  6 

T.  R.  Marrero  and  E.  A.  Mason,  "Gaseous  Diffusion  Coefficients",  J.  Phys. 
Chem.  Ref.  Data  1,  3-118,  (1972). 
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TABLE  2.  THERMODYNAMIC  COEFFICIENTS 
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TABLE  3.  TRANSPORT  COEFFICIENTS 
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To  complete  the  set  of  input  coefficients  expressions  for  A3,  pD23 
and  pDi3  must  be  found.  To  accomplish  this  we  first  obtain  from  vis¬ 
cosity  data,  Lennard -Jones  (L-J)  parameters  for  0  and  02  and  Stockmayer 
parameters  for  O3.  Then  theoretical  formalisms8 > 2?  are  USed  to  deter¬ 
mine  values  for  each  transport  coefficient  as  a  function  of  temperature. 
Finally  a  fit  of  the  form  AT8  is  made  to  the  values  and  listed  in  Table 
3.  An  outline  of  the  details  follows. 

The  L-J  parameter  ei/k  is  found  from  the  ratio  of  the  values  of  the 
viscosity  at  two  temperatures.28  The  atomic  oxygen  viscosity  is  taken 
fgon^the  expression  of  Dalgarno  and  Smith,24  i.e.,  ni  =  3.34  x  10-6 
T  ’  gm-cm-l-s-4,  100K  <  T  <  2000K.  The  molecular  oxygen  viscosity  is 
taken  from  the  tabular  data  of  Hanley  and  Ely.25  (We  elect  to  determine 
our  own  set  of  L-J  parameters  for  02  rather  than  use  Svehla 's29  because 
Hanley  and  Ely  noted  internal  inconsistencies  in  certain  viscosity 
measurements  that  were  available  to  Svehla) .  Having  found  ei/k  we  can 
determine  the  L-J  parameter  from  a  direct  substitution  into 

266.93  x  10  (M.T)  '  f  (T.) 

O-(T)  =  - ~ — j- i - * - —  gm-cm  -s”»  (11) 

V  ft  *  (T.) 

In  practice  we  fix  300K  as  the  first  temperature  and  compute  sets  of  L-J 
parameters  for  each  second  temperature  in  the  interval  400K-2000K  in 
steps  of  100K.  The  L-J  parameters  that  give  the  best  fit  to  the  viscosity 
data  in  a  least  squares  sense  are  chosen  and  listed  in  Table  4.  Over  the 
range  300K-2000K  our  L-J  parameters  for  molecular  oxygen  reproduce  the 
tabular  viscosity  of  Hanley  and  Ely  to  within  2%  while  Svehla' s  L-J 
parameters  (02/^  =  106. 7K  and  a~  =  3.467A)  reproduce  the  same  data  within 
4%.  z 


The  determination  of  the  Stockmayer  parameters  for  ozone  is  compli¬ 
cated  by  the  fact  that  we  have  available  only  the  viscosity  of  ozone- 
molecular  oxygen  mixtures.  Tretter5^  records  five  measurements  of  the 
mixture  viscosity  of  nearly  equal  molar  volumes  of  O2  and  O3  at  194K, 


27 

R.  C.  Reid  and  T.  K.  Sherwood 3  The  Properties  of  Gases  and  LiguidSj 
2nd  Edition 3  McGraw-Hill3  NY3  (1966). 

28 

J.  0.  Hirschfelder3  C.  F.  Curtis  and  R.  B.  Bird 3  op,  citj_3  p.  562. 

29 

R.  A.  Svehla3  " Estimated  Viscosities  and  Thermal  Conductivities  of 
Gases  at  High  Temperatures" ,  NASA  Tech.  Report  R-132 3  Lewis  Res.  Ct.3 
Cleveland ,  Ohio3  (1962). 

30 

V.  J.  Tretter3  An  Ex^ervnental  Determination  of  The  Viscositii  of  Gaseous 
Ozonej  M.  S.  thesis 3  U.  of  Idaho 3  (1966) 3  unpublished. 
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273K  and  295K.  The  procedure  we  have  used  to  recover  the  ozone  Stockmayer 
parameters  from  Tretterfs  data  is  as  follows. 

The  expressions  for  the  mixture  viscosity  as  given  by^l  are  functions 
of  the  L-J  (Stockmayer)  parameters,  the  temperature,  molar  composition 
and  certain  combining  rules.  For  the  O2/O3  mixture  the  L-J  parameters 
for  molecular  oxygen  are  determined  above  and  the  temperature  and  com¬ 
position  are  fixed  by  the  experimental  conditions.  The  empirical  polar/ 
non-polar  combining  rules^  are 


°23  “  0-5  *-ct2  +  a3-)  ^  ^  (12) 


and 


'23 


-  r  ,1/2  (.2 

-  (e2e3)  5  : 


where  5  =  1  +  0.25  a2  y32  [a23  c?33  (£2^3) 1/2] _1  •  (The  value  used  for 
Th®  5)ean  dipole  polarizability  of  molecular  oxygen,  a2,  is  1.6  x  lO-*"4 
cm  .  >  ^  The  ozone  dipole  moment  is  taken  as  the  average  of  the  four 

values  listed  by  McCellan;34  i.e.,  p3  =  0.57  x  10-18  dyne  l/2_cm2.)  The 
only  unknowns  in  the  mixture  viscosity  expressions  are  the  Stockmayer 
parameters  for  ozone.  A  search  is  made  over  the  two  dimensional  space 
of  the  Stockmayer  parameters,  pairs  of  parameters  are  selected  and  a 
mixture  viscosity  computed.  A  computed  value  is  compared  with  the  cor¬ 
responding  measured  value  and  a  root-mean-square  (RMS)  error  is  deter¬ 
mined  for  the  five  temperature/composition  conditions  of  the  measurements. 
The  minimum  RMS  error  then  determines  the  ozone  Stockmayer  parameters 
listed  in  Table  4.  This  procedure  with  Svehla's  L-J  parameters  for  mole¬ 
cular  oxygen  results  in  no  essential  change  in  the  Stockmayer  parameters 
for  ozone. 

The  Stockmayer  parameters  for  ozone  are  used  in  the  appropriate 
form  of  equation  (11)  to  find  the  ozone  viscosity  as  a  function  of  tem¬ 
perature.  We  find  P3(T)  =  4.24  x  10"6t0.65  gm-Cm-l-s-l.  Finally  A3  is 
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J.  0.  Rirschf elder,  C.  F.  Curtis  and  R.  B.  Bird,  op.  git. ,  p.  530. 

32 

J.  0.  Rirsehf elder,  C.  F.  Curtis  and  R.  B.  Bird,  git. ,  p.  600. 

33 

J.  0.  Rirschf 'elder ,  C.  F.  Curtis  and  R.  B.  Bird,  op.  cit^ ,  p.  950. 

S4A.  L.  McCellan,  Tables  of  Ex£erimental  Dipole  Moments ,  W.  H.  Freeman 
and  Co.,  San  Francisco,  (1963). 
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TABLE  4.  LENNARD-JONES/STOCKMAYER  PARAMETERS 


O 


lecies 

a  (A) 

e/k  (K) 

0 

2.947 

127.2 

°2 

3.372 

128.7 

°3 

4.69 

39.4 
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found  from  t)3  by  means  of  Reid  and  Sherwood*  s27  recommended  relationship 
between  heat  conductivity  and  viscosity,  the  Bromley  analysis. 

Using  the  parameters  listed  in  Table  4  and  the  empirical  combining 
laws  given  by  equation  (12)  the  diffusion  coefficients  pDi3  and  pD23 
are  computed  as  functions  of  temperature  from  the  expression^ 3 27 


2.628  x  10 *  3 *  [T3 
2 

a .  . 
13 


(M.  +  M.)/  2  M.M.]1//2 
1  3  1  3 

11  * 

n  ’  (T.  .) 

13 


«.  2  -1 
atmos-cm  -s 


(13) 


Measurements  for  comparison  are  all  but  non-existent.  Hildebrand35 
has  discovered  an  empirical  correlation  that  links  the  heat  capacity  of 
a  gas  with  its  viscosity  at  298K.  Using  the  ozone  heat  capacity  of  9.4 
cal-mole-l-K-1  at  298K22  and  Hildebrand's  correlation,  we  find  ri 3 C 2 98 K) 

=  164  x  10-6  gm-cm-ls-l.  The  value  computed  from  our  expression  for 
03  (298K)  is  168  x  10_t)  gm-cm--*- -s“l .  Yastrebov3^  has  estimated  A3  at 
room  temperature  to  be  about  three-quarters  that  of  A2.  From  Table  3 
we  find  0.75  X2  (300K)  =  4.82  x  10~5  cal-cm-l-s-i-K-l  and  A3  (300K) 

=  4.75  x  10"5  cal-cm_l-s-l_K-l .  Houzelot  and  Villermaux  are  reported"^ 
to  have  measured  U23  =  0.15  cm2-s_1  at  one  atmosphere  and  25°C.  From 
Table  3  we  find  D23(298K)  =  0.155  cm2-s_1. 

Another  independent  check  can  be  made.  From  the  L-J  parameters 
for  0  and  O2,  (pDi2)gj  can  be  formed  using  the  appropriate  analogue  of 
equation  (13)  and  equation  (12)  (yi  =  y2  =  0  implies  £  =  1).  This 
computation  yields  the  fit  (pUi2)Lj  =  2.01  x  10”5  Tl-687  atmos-cm2-s _1 
over  the  range  300K  -  2000K.  The  ratio  (pDi2)gj  to  pDi2  shows  a 
monotonic  decrease  from  0.92  at  300K  to  0.78  at'  2000K.  This  comparison 
provides  us  with  a  measure  of  the  efficacy  of  the  procedures  outlined 
above . 


?°J.  H.  Hildebrand >  Mole.  Phys.  35 ^  519-523,  (1978). 

3  6 

V.  V.  Yastrebov,  "The  Physical  Chemistry  of  Concentrated  Ozone,  VIII. 
Thermal  Propagation  of  Flame  in  Gaseous  Ozone  Mixtures" ,  Russian  J.  of 
Phys.  Chem .,  34j,  21-23,  (1960). 

37 

J.  L.  Houzelot  and  J.  Villermaux,  "Measurement  of  the  Molecular 
Diffusitivity  of  Ozone  Oxygen",  C.  R.  Acad.  Sci.  Ser.  C.  1971,  273, 
258-60.  As  quoted  by  Chem.  Abstracts  74,  252,  (1971). 
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V 


Equation  (3)  requires  an  expression  for  the  heat  conductivity  of 
the  mixture,  X.  We  have  used3** 


N  N  -1 

X  =  0.5  E  X.X.  +  E  X./X. 

.  .  li  .  .  l  i 


V.  RESULTS  AND  DISCUSSION 

Figures  1-4  show  our  computed  profiles  for  four  initial  ozone  mole 
fractions  (IOMF's).  For  each  IOMF  Tu  =  300K  and  the  total  pressure  is 
one  atmosphere.  We  do  not  know  of  any  experimental  profiles  for  com¬ 
parison  and  so  by  our  definition  (see  Introduction)  the  validation  of 
the  ozone  flame  remains  to  be  done.  However,  experimental  burning  velo¬ 
cities  are  available. 

Figure  5  shows  the  burning  velocity  as  a  function  of  IOMF.  Our 
computed  values  of  the  velocities  are  plotted  as  circles  and  the  experi¬ 
mental  data  of  Streng  and  Grosse39  are  shown  as  squares.  The  solid  line 
is  their  expression  for  the  burning  velocity  and  is  given  by 


Su  =  563  Xq  -  88.8  cm/s.  (14) 


Table  5  lists  our  computed  and  the  experimental  values  for  the  burning 
velocity.  This  right-hand,  column  shows  that  between  IOMF  of  0.25  and 
1.00,  our  computed  values  are  5  to  28%  greater  than  the  experimental 
values.  Though  no  exact  figures  are  given,  Streng  and  Grosse  state  that 
for  measurements  less  than  0.3  IM0F  their  experimental  results  are  less 
accurate.  Thus  the  38%  discrepancy  for  0.20  IOMF  is  not  necessarily  a 
defect  of  the  model. 

Tests  were  made-*-^  at  IOMF  =  0.25  and  1.0  to  examine  the  sensitivity 
of  the  computed  burning  velocity  to  changes  in  the  number  of  breakpoints, 
in  the  polynomial  order  of  the  spline  functions,  KORD,  and  in  the  number 
of  continuity  conditions  imposed  at  the  breakpoints,  NCC.  For  each  case 
we  started  the  model  from  a  standard  profile  and  ended  the  exercise  at 
a  time  sufficient  to  ensure  that  steady  state  had  been  achieved.  For 
our  standard  conditions  KORD  =  6  and  NCC  =  5  a  doubling  of  the  number  of 


38 

Y.  S.  Touloukian,  P.  E.  Liley  and  S.  C.  Saxena,  Thermoph^sieal  Properties 
of  Matter^  Vol.  3,  Thermal  Conductivity,  Nonmetallic  Liquids  and  Gases), 
IF I /Plenum,  NY-Washington,  (1970),  p  45a. 

39 

A.  G.  Streng  and  A.  V.  Grosse,  "The  Ozone  to  Oxygen  Flame",  Sixth 
Symposium  (International)  on  Combustion ,  Reinhold  Publishing  Company, 
(1957),  pp  264-273. 


25 


1400 


(»)  3amva3dW3i 


O  o  O  O  O  O  o 

O  o  o  o  o  o 

CM  O  °0  O  CN 


26 


-1.0  -0.5  0 

DISTANCE  (mm) 

Figure  1.  Temperature  and  Species  Profiles  for  IOMF  =  0.25,  Tu  =  300K  and  Total  Pressure 
1  atmosphere. 
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Figure  2.  Temperature  and  Species  Profiles  for  IOMF  =  0.50,  Tu  =  300K  and  Total  Pressure 
1  atmosphere. 
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Figure  3.  Temperature  and  Species  Profiles  for  IOMF  =  0.75,  Tu  =  300K  and  Total  Pressure 
1  atmosphere. 
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Figure  4.  Temperature  and  Species  Profiles  for  IOMF  =  1.0,  Tu  =  300K  and  Total  Pressure 
1  atmosphere. 
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INITIAL  OZONE  MOLE  FRACTION 

Figure  5.  Burning  Velocity  vs.  I0MF.  Squares  are  the  experimental  points.  Circles  are  computed 
from  our  model.  The  line  is  Streng  and  Grosse's  expression  for  the  burning  velocity. 


TABLE  5.  COMPARISON  OF  OUR  COMPUTED  BURNING  VELOCITIES 
WITH  THE  MEASUREMENTS  OF  STRENG  AND  GROSSE 


Initial  Ozone 

Mole  Fraction 

Burning  velocity  (cm/s) 

This  work  Streng  $  Grosse* 

Ratio 

1.00 

497 

474 

1.05 

0.75 

396 

333 

1.19 

0.50 

248 

193 

1.28 

0.25 

64 

52 

1.23 

0.20 

33 

24 

1.38 

* Determined,  from  equation  (14). 
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breakpoints  produced  a  maximum  difference  of  0.3%  in  the  burning  velocity 
computed  using  a  given  species  (see  equation  (10)).  By  way  of  comparison 
the  maximum  species  to  species  difference  in  the  burning  velocity  was 
0.4%.  For  IOMF  =  1.0  we  held  the  number  of  breakpoints  fixed  at  70  and 
subject  to  the  condition  NCC  <  KORD,  varied  both  KORD  and  NCC  over  the 
range  3  <_  KORD  8  and  2  -c  NCC  7.  The  maximum  change  in  the  burning 
velocity  computed  from  the  atomic  oxygen  (02,03)  profile  was  2.6%  (1%). 

A  similar  test  for  IOMF  =  0.25  showed  a  maximum  change  of  1.9%  (0.3%). 

We  infer  from  the  observed  insensitivity  that  our  standard  conditions 
are  sufficient  to  produce  a  numerically  reliable  solution.  We  also  found 
that  our  computational  efficiency  could  have  been  improved  over  our 
standard  conditions  by  reducing  both  KORD  and  NCC  by  one  or  two. 

Table  6  shows  the  percentage  change  in  velocities  of  IOMF's  of  0.25 
and  1.00  for  a  given  change  in  each  of  the  kinetic  and  transport  coeffi¬ 
cients.  (The  kinetic  coefficients  have  been  altered  in  pairs  so  that 
the  equilibrium  constant  remains  the  same).  The  factor  of  1.2  is 
arbitrary;  but  we  feel  representative  for  the  input  coefficients  except 
for  k3  (factor  of  two)  and  for  A  2  (<_  6%) .  Because  of  non-linear  effects 
we  do  not  expect  these  percentage  changes  to  be  quantitatively  scaleable 
to  other  changes  in  the  input  coefficients;  nonetheless,  we  shall  use 
them  to  indicate  trends  in  the  relative  importance  of  each  coefficient. 

Warnatz-^  has  recently  and  independently  assembled  a  finite  differ¬ 
ence  code  that,  like  ours,  requires  species  dependent  input  coefficients. 
He  has ^ found  that  the  discrepancies  in  the  model  results  of  Dixon-Lewis , 1 
Wilde, ^  Bledjian^  and  Campbell^  are  traceable  to  the  different  expressions 
for  the  input  coefficients  each  used.  In  the  light  of  this  discussion,^ 
we  confine  our  remarks  to  a  comparison  of  our  model  results  with  those 
of  Warnatz. 


Table  7  shows  the  comparison  of  the  burning  velocity  for  five  IOMF's. 
Warnatz'  values  were  read  from  his  Figure  6,  curve  labelled  "M".  The 
right  hand  column  shows  the  ratio  of  Warnatz'  computed  burning  velocity 
to  ours.  It  is  seen  that  the  two  model  results  agree  within  ±  12%  over 
the  range  0.20  <_  IOMF  <_  1.00.  This  agreement  does  not  imply  that  the 
sets  of  input  coefficients  used  in  the  respective  models  are  equivalent. 

Figure  6  plots  the  ratio  of  the  values  of  Warnatz'  input  coefficients 
to  our  corresponding  values  over  the  temperature  range  1000K  ■<  T  <  2400K. 
The  striking  features  of  this  plot  are  the  ratios  for  ki  and  k2-  _The 
value  for  the  kj  ratio  falls  from  0.92  at  1000K  to  0.39  at  2400K  while 
k2  rises  from  1.46  to  1.76  over  the  same  temperature  range.  Since  we 
had  both  used  the  same  source^  for  k3,  its  ratio  is  unity  and  is  not 
shown  in  the  Figure.  (The  reverse  of  all  three  reactions  plays  no  sig¬ 
nificant  part  in  the  determination  of  the  burning  velocity  and  they  are 
not  further  discussed.) 

The  ratios  of  the  values  of  the  transport  coefficients  are  seen  to 
be  nearly  independent  of  temperature  (both  Warnatz  and  we  have  used 
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TABLE  6.  PERCENTAGE  CHANGE  IN  THE  COMPUTED  BURNING  VELOCITY 
AT  INITIAL  OZONE  MOLE  FRACTIONS  OF  0.25  AND  1.00  WHEN  EACH 
INPUT  COEFFICIENT  IS  MULTIPLIED  BY  A  FACTOR  OF  1.2 


Input 

Coefficient 

%  Change  in  Burning  Velocity 

0.25  1.00 

kl 

8.6 

4.6 

k2 

0.8 

4.4 

k3 

0.0 

0.0 

X1 

0.2 

0.4 

X2 

11.8 

3.6 

X3 

0.9 

2.5 

pD12 

2.7 

1.9 

pD13 

0.5 

4.5 

pD23 

-5.7 

-3.6 

TABLE 

7.  COMPARISON  OF  WARNATZ'  COMPUTED 
BURNING  VELOCITIES  WITH  OURS 

Initial  Ozone 

Mole  Fraction 

Burning  Velocity  (cm/s) 

Warnatz  This  work 

Ratio 

1.00 

445 

497 

0.90 

0.75 

350 

396 

0.88 

0.50 

225 

248 

0.91 

0.25 

65 

64 

1.02 

0.20 

37 

33 

1.12 

33 


RATIO  INPUT  COEFFICIENTS 


Figure  6.  Ratio  of  the  Values  of  Warnatz1  and  Our  Input  Coefficients; 

e.g.,  the  ratio  k]  .Warnatz/M ,thi s  paper  is  given  by  the 
curve  identified  by  k]  . 
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similar  functional  forms)  and  the  respective  values  of  each  ratio  with¬ 
in  ±  20%  of  each  other.*  The  important  point  is  that  the  differences 
in  transport  coefficients  are  much  less  than  those  associated  with  the 
kinetic  coefficients  ki  and  k2. 

If  the  kinetic  coefficients  for  k^  and  k2  differ  by  so  much  how 
can  the  burning  velocities  listed  in  Table  7  be  in  essential  agreement? 

The  answer  is  suggested  by  the  sensitivities  of  the  burning  velocity  to 
changes  in  ki  and  k2  as  listed  in  Table  6. 

For  the  1.00  IOMF  the  influence  ki  and  k2  upon  the  burning  velocity 
are  about  equal;  and,  in  going  from  one  model  to  the  other,  the  high 
temperature  values  of  ki  and  k2  change  in  opposite  directions  (see 
Figure  6) .  Thus  in  this  case  there  is  a  tendency  for  the  individual 
effects  of  ki  and  k2  on  the  burning  velocity  to  offset  each  other. 

This  tendency  can  be  checked  in  the  following  manner.  Using  our 
model  results  for  the  burning  velocity  at  an  IOMF  of  1.00  as  a  fiducial, 
we  successively  substitute  into  our  model  Warnatz'  expressions  for  k]_, 
then  for  k2,  and  finally  for  both  ki  and  k2  and  compute  the  burning 
velocity  in  each  case.  As  can  be  seen  in  Table  8  the  burning  velocity 
so  computed,  first  goes  down  relative  to  the  fiducial,  then  up  and  finally 
shows  the  mutually  offsetting  effects  of  these  changes  in  k]^  and  k2-  As 
an  added  bonus  the  burning  velocity  of  459  cm/s  obtained  by  substituting 
Warnatz'  kinetic  coefficients  for  k]^  and  k2  into  our  code  (see  Table  8) 
can  be  compared  to  Warnatz'  own  computed  value  of  445  cm/s  to  provide  a 
measure  of  the  collective  effects  of  differing  transport. 

Table  6  shows  that  for  the  0.25  IOMF  case,  the  burning  velocity's 
sensitivity  to  changes  in  kinetic  coefficients  is  dominated  by  ki. 

Figure  6  shows  that  Warnatz'  and  our  values  for  ki  are  much  closer  at 
the  lower  temperatures. 

In  comparing  profiles  there  are  sensible  differences  in  some  model 
results.  For  1.00  IOMF,  Figures  7  and  8  show  a  comparison  of  the  atomic 
oxygen  and  temperature  profiles,  respectively.  (Warnatz'  values  were 
read  from  his  Figure  7).  For  ease  in  viewing,  the  curves  have  been 
arbitrarily  displaced  from  each  other  along  the  distance  axis  in  both 
figures.  Figure  7  shows  that  within  the  burned  region,  behind  the  flame- 
front,  i.e.,  -0.05  to  -0.10mm,  our  computed  value  for  the  atomic  oxygen 
mole  fraction  lies  about  a  factor  of  two  greater  than  Warnatz'.  Figure 
8  shows  that  over  this  same  region  our  computed  temperature  lies  about 
350K  below  Warnatz' . 


We  note  that  the  L-J  and  Stockmayer  parameters  used  by  Warnatz  are  quite 
different  than  those  used  here.  The  fact  that  the  transport  values  oan 
be  represented  by  different  pairs  of  parameters  has  been  discussed.  2, 27 
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TABLE  8.  SUCCESSIVE  SUBSTITUTION  OF  WARNATZ'  (W)  KINETIC 
COEFFICIENTS  FOR  OURS  (0)  AT  AN  INITIAL  OZONE  FRACTION  OF  UNITY 


k 


1 


Burning  Velocity 


0  0  497  cm/s 

W  0  427 


0 


W 


S53 


W 


w 


459 
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Figure  7.  Our  and  Warnatz'  Calculated  Atomic  Oxygen  Profiles  for  IOMF  =  1.0.  (These  curves  are 

arbitrarily  displaced  for  ease  in  viewing.)  Dashed  profile  results  for  the  substitution 
of  Warnatz1  expressions  for  k,  and  k9  into  our  model. 
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Figure  8.  Our  and  Warnatz'  Calculated  Temperature  Profiles  for  IOMF  =  1.0.  (These  curves  are 

arbitrarily  displaced  for  ease  in  viewing.)  Dashed  profile  results  from  the  substitution 
of  Warnatz1  expressions  for  k,  and  k0  into  our  model. 


We  now  consider  the  causes  of  the  profile  differences.  In  the 
burned  region,  where  the  ozone  concentration  and  the  species  and  thermal 
gradients  are  small,  the  only  reaction  of  importance  is  the  atomic 
oxygen  recombination  reaction.  (We  find  the  effect  of  k  completely 
negligible  which  is  equivalent  to  taking  k_3  =  0  as  was  done  by  Warnatz.) 

As  noted  previously,  both  models  employ  the  same  expression^ 7  fQr  k3. 

Thus,  the  causes  of  these  differences  in  the  burned  region  must  be 
sought  in  the  flame-front.  In  this  region,  as  we  have  seen,  the  main 
differences  in  the  model  input  coefficients  are  the  expressions  for  k^ 
and  k2.  Substitution  of  Warnatz'  expressions  for  k^  and  k2  into  our  code 
yields  the  dashed  line  in  both  Figures  7  and  8.  This  substitution  shows 
that  the  preponderance  of  the  differences  in  these  profiles  is  due  to 
the  different  expressions  for  k]_  and  k2-  In  other  words  these  different 
expressions  lead  to  different  atomic  oxygen  superequilibrium  concentrations 
and  the  relaxation  to  the  final  equilibrium  state,  via  reaction  (III), 
is  so  slow  that  these  differences  are  effectively  frozen. 


The  fact  that  both  models  employ  the  same  expression  for  k3  permits 
us  to  check  on  the  models'  mutual  consistency  in  the  burned  region.  To 
wit,  the  atomic  oxygen  difference  between  model  results  as  seen  in 
Figure  is  sufficient  to  account  for  the  temperature  difference  between 
the  models'  results  as  seen  in  Figure  8. 


For  0.20  IOMF ,  the  values  of  the  controlling  input  coefficients  for 
the  two  models  are  similar.  We  expect  and  find  similar  species  and 
temperature  profiles  (not  shown). 

We  have  found  that  the  area  of  a  greatest  discrepancy  in  the  input 
values  used  for  these  two  models  lies  in  the  expressions  used  for  the 
two  rate  coefficients,  k^  and  k2*  We  have  recently  critically  evaluated 
the  available  high  temperature  experimental  data  for  the  ozone  decompo¬ 
sition  reaction  and  the  expression  used  herein  is  consistent  with  all 
the  direct  experimental  data  known  to  the  authors. 18  (jn  order  to  use 
the  various  experimental  data  to  obtain  the  ozone  decomposition  expression 
we  have  assumed  that  the  relative  third  body  efficiencies  are  independent 
of  temperature) . 

The  expression  used  for  k2  is  another  story.  Warnatz  has  developed 
and  used  his  own  expression^  and  we  have  employed  Hampson ' s . 21  To  use 
them  in  our  codes  we  have  both  assumed  that  the  respective  expressions 
are  valid  beyond  the  upper  limit  of  applicability  of  each,  i.e.,  1000K. 

In  order  to  distinguish  which  expression  for  k2,  if  either,  is 
correct,  measurements  and/or  ab  initio  calculations  of  the  rate  coeffi¬ 
cient  for  reaction  (II)  are  required  at  high  temperatures,  say  in  the 
range  1500-2000K.  Alternately,  the  computed  differences  in  the  values 
for  atomic  oxygen  and  the  temperature  in  the  burned  region  at  an  IOMF 


J.  Warnatz 3  "Berechnung  der  Flammengeschwindigkeit  und  der  Struktur 
von  laminaren  flachen  Flarrmen ",  Habilitationssohrift 3  Darmstadt ,  1977. 
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of  1.00  appear  to  be  large  enough  that  profile  measurements  above  a 
flame  may  be  sufficient  to  distinguish  between  the  two  expressions. 


VI.  SUMMARY 

We  have  shown  that  this  model  and  its  input  parameters  (see  Tables 
1,  2,  and  3)  predict  burning  velocities  that  are  in  reasonable  agreement 
with  the  measurements  of  Streng  and  Grosse39  and  with  the  computations 
of  Warnatz.3  We  have  also  demonstrated  that  agreement  with  burning 
velocities,  even  over  a  wide  range  of  initial  ozone  mole  fractions,  is  a 
necessary  but  not  sufficient  condition  to  ensure  that  the  input  coeffi¬ 
cients  are  realistic.  We  have  shown  that  profile  measurements  are  vital 
to  test  the  model's  input  coefficients  and  so  validate  the  model. 
Finally,  by  a  comparison  of  computed  profiles  we  have  indicated  the  need 
to  measure  or  to  calculate  high  temperature  values  for  k2. 
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GLOSSARY 


Pi 


h. 

1 


k. 

i 


N 

M. 

i 

m 

c 

P 

R 


R. 

l 


* 

T. 

l 


V. 

l 


X. 

1 


=  Heat  capacity  at  constant  pressure,  cal-mole-1-^1 
=  s  CpiYi»  the  specific  heat  of  the  mixture,  cal-gm"1-K”1 
=  specific  heat  of  ith  species,  cal-gm'1-K-1 

—  1 

=  binary  diffusion  coefficient  for  species  k  and  j,  cm  -s 
=  enthalphy,  cal-mole”'*' 

=  specific  enthalpy,  cal-gm”'*' 

=  Boltzmann's  constant,  =  1.3806  erg  -  K-1 

=  rate  coefficient  for  the  ith  reaction,  in  centimeter-mole-seconds 
units 

=  number  of  species 

=  molecular  weight  of  ith  species,  gm-mole"1 

-2  -1 

=  mass  flux,  gm-cm  -s 
=  total  pressure,  atmosphere 

=  Gas  constant,  =  1.9872  cal-mole" 1K_1 ,  =  82.05  cm3-atmos-mole"1-K_1 

=  Rate  of  production  of  ith  species  by  chemical  reactions, 
mole-cm-3-s-l 

=  burning  velocity,  cm-s-1 
=  temperature,  K 
=  TCc./k)-1 

=  temperature  of  unburned  mixture,  K 
=  temporal  coordinate,  s 
=  fluid  velocity,  cm-s"^ 

=  diffusion  velocity  of  ith  species,  cm-s"1 
=  spatial  coordinate,  cm 
=  mole  fraction  of  ith  species 
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Yi  =  mass  fraction  of  ith  species 

Yiu  =  mass  fraction  of  ith  species  in  the  unburned  mixture 

ej/k  =  L-J  or  Stockmayer  parameter,  K 

**  viscosity  of  ith  species,  gm-s  ^-cm  ^ 

A  =  mixture  conductivity,  cal-cm"1-s~1-K_1 

Xi  =  conductivity  of  ith  species,  cal-cm"1-s"1-K"1 

Q  O 

°i  =  L-J  or  Stockmayer  parameter,  10  cm  (=  A) 

p  =  fluid  density,  gm-cm"3 

^  -  Lagrangian  coordinate  (see  equation  6) ,  gm-cm”"* 

01,2  2,2  _  ... 

5  ~  collision  intergrals  (tabulated  in  reference  8) 
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